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ABSTRACT 

The Visible Infrared Imager Radiometer Suite (VIIRS) instrument onboard the Suomi National Polar-orbiting 
Partnership (SNPP) satellite was launched on 28 October 2011. The VIIRS has 5 imagery spectral bands (I-bands), 16 
moderate resolution spectral bands (M-bands) and a panchromatic day/night band (DNB). Performance of the VIIRS 
spatial response and band-to-band co-registration (BBR) was measured through intensive pre-launch tests. These 
measurements were made in the non-aggregated zones near the start (or end) of scan for the I-bands and M-bands and 
for a limited number of aggregation modes for the DNB in order to test requirement compliance. This paper presents 
results based on a recently re -processed pre-launch test data. Sensor (detector) spatial impulse responses in the scan 
direction are parameterized in terms of ground dynamic field of view (GDFOV), horizontal spatial resolution (HSR), 
modulation transfer function (MTF), ensquared energy (EE) and integrated out-of-pixel (IOOP) spatial response. 
Results are presented for the non-aggregation, 2-sample and 3 -sample aggregation zones for the I-bands and M-bands, 
and for a limited number of aggregation modes for the DNB. On-orbit GDFOVs measured for the 5 I-bands in the scan 
direction using a straight bridge are also presented. Band-to-band co-registration (BBR) is quantified using the pre- 
launch measured band-to-band offsets. These offsets may be expressed as fractions of horizontal sampling intervals 
(HSIs), detector spatial response parameters GDFOV or HSR. BBR bases on HSIs in the non-aggregation, 2-sample 
and 3 -sample aggregation zones are presented. BBR matrices based on scan direction GDFOV and HSR are compared 
to the BBR matrix based on HSI in the non-aggregation zone. We demonstrate that BBR based on GDFOV is a better 
representation of footprint overlap and so this definition should be used in BBR requirement specifications. We propose 
that HSR not be used as the primary image quality indicator, since we show that it is neither an adequate representation 
of the size of sensor spatial response nor an adequate measure of imaging quality. 

Keywords: Suomi NPP, VIIRS, band-to-band co-registration, ground dynamic field of view, spatial response, 
ensquared energy, integrated out-of-pixel spatial response, day/night band (DNB). 


1. INTRODUCTION 

The NASA/NOAA Visible Infrared Imager Radiometer Suite (VIIRS) instrument onboard the Suomi National 
Polar-orbiting Partnership (SNPP) satellite was launched on 28 October 2011. SNPP was put under management of the 
NASA acquired and NOAA managed Joint Polar Satellite System (JPSS) after restructuring of the National Polar- 
orbiting Environmental Satellite System (NPOESS) and it was renamed from the NPOESS Preparatory Project (NPP) 
after its successful launch and check out in early 2012 [11 . Copies of VIIRS instrument will be on the follow-on NOAA 
Joint Polar Satellite System (JPSS) in an afternoon orbit (1:30 PM local time of the ascending node crossing the 
equator). As a payload originally developed for the NPOESS, the VIIRS instrument inherits the triple-purpose (civilian, 
military and science applications) design 12 ’ that balances requirements such as fine imagery resolution, low noise and 
broad spectral coverage. 
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VIIRS has 5 imagery resolution bands (bands II to 15) with 32 detectors each, 16 moderate resolution bands (bands Ml 
to M16) and a panchromatic day-night band (DNB) with 16 detectors each. The spectral band M16 has two physical 
bands M16A and M16B that are time delay integrated spatially but tested individually pre-launch on the ground. DNB 
has three gain stages in four charge-coupled device (CCD) arrays: high gain A-side (HGA), high gain B-side (HGB), 
mid-gain stage (MGS) and low-gain stage (LGS). These gain stages are selected in on-orbit operations but tested 
individually pre-launch. These bands collect energy from the earth in a cross-track (scan) direction through a rotating 
telescope assembly (RTA) and a half angle mirror (HAM) that effectively de-rotates the rays from the RTA 1 ' 1 . As the 
satellite circles the earth in the near -polar orbit and the RTA scans across the track, VIIRS sweep a swath of ~3,000 km 
of the earth in a field of regard of +/- 56°. The M-bands have a nominal horizontal sampling interval (HSI) of 750 m at 
nadir, the I-bands have a nominal HSI of 375 m at nadir, and the day-night band (DNB) has a near-constant nominal 
750 m HSI throughout the scan. Due to the effects of slant distance and earth curvature, the HSI from each detector 
grows from nadir to the edges of scan and ground coverage overlaps more and more towards the edge of scan (“bow- 
tie” effect). An along-scan aggregation scheme is employed to limit the growth of pixels in the scan direction. In the 
nadir region of the scan, three samples are aggregated (called 3x1 aggregation or Agg3xl). No aggregation (called lxl 
aggregation or Agglxl) is performed near the edges of the scan and two samples are aggregated (called 2x1 
aggregation or Agg2xl) in the middle of the scan. Also, a bow-tie deletion scheme is applied to trim readings of one 
(two) and two (four) M (I) band detectors at the edges of every band detector array in the Agg2xl and Agglxl zones, 
respectively. For the DNB, 32 distinctive sample aggregation and bow-tie “compression” (contiguous CCD detectors 
selected to form 16 “super” detectors close to the center in the track direction) is used to keep the ground sampling 
interval near constant throughout the scan with Aggregation Mode 1 at nadir and Aggregation Mode 32 at the edges of 
scan 11 ’ 6 ' 7 '. Performance of spatial response and band-to-band co-registration (BBR) was measured through intensive pre- 
launch tests and the results have been presented in terms of requirement compliance 14 ’ 5 ' 81 . 

This paper presents results after re -processing the pre-launch data from thermal vacuum tests at the nominal 
temperature performance plateau. These ground tests were set in the diagnostic mode, meaning non-aggregation for all 
detectors for the I-bands and M-bands. Line spread functions (LSFs) and band-to-band offsets in the Agg2xl and 
Agg3xl zones are constructed analytically from the non-aggregation test data. For the DNB spatial responses, 
aggregation modes 18, 22 and 31 were tested for all four CCD arrays. Section 2 presents the along scan LSF parameters 
in terms of dynamic field of view (DFOV), horizontal spatial resolution (HSR), modulation transfer function (MTF), 
ensquared energy (EE) and integrated out-of-pixel (IOOP) spatial responses. Results are presented for the three zones 
for the I-bands and M-bands, and for the three aggregation modes for the DNB. On-orbit GDFOVs for the 5 I-bands in 
the scan direction using a straight bridge are also presented in this Section. Section 3 presents the I-band and M-band 
band-to-band co-registration (BBR) matrices based on HSIs in the three aggregation zones. BBR matrices based on 
scan direction GDFOV and HSR are compared to the BBR matrix based on HSI in the non-aggregation zone in this 
Section as well. Concluding remarks are presented in the last Section. 


2. VIIRS SENSOR SPATIAL RESPONSES 

The VIIRS sensor spatial impulse responses were measured for every detector before launch. Since the detectors were 
fabricated in an optically apparent rectangular shape, their full impulse responses were determined by combining two 
orthogonal line spread functions (LSFs) measurements. These LSFs were measured during the thermal vacuum (TVac) 
testing. The tests included LSFs measurements in the scan and track directions at the cold and nominal performance 
temperature plateaus, and some special tests at the hot performance plateau. The scan direction LSFs presented here are 
from the test at the nominal performance temperature plateau, representing the nominal operating conditions in space. 
The LSFs in the track direction are mostly squares 14 ’ 81 , with their width nearly the same as the horizontal sampling 
interval (HSI). 

The ground tests were performed in the diagnostic mode, meaning none of the detectors were aggregated. Band Ml 6 is 
time delay integrated from bands M16A and M16B, but only band M16B was measured and included in this paper. The 
LSFs in the Agg2xl and Agg3xl zones are derived from the LSFs in the Agglxl zone. Angular sampling intervals 
(ASIs) in the scan direction are constants for I-bands and M-bands, respectively at 156 and 312 pradians in the Agglxl 
zone. The ASIs in the Agg2xl and Agg3xl zones are respectively two and three times as large as those in the Agglxl 
zone. In the track direction, ASIs are constants throughout the scan, at 445.5 and 891 Liradians for I-bands and M-bands, 
respectively 141 . For DNB, the track and scan ASIs varies with scan angle in 32 aggregation modes, from 894 uradians 
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for the aggregation mode 1 at nadir to 149 pradians for the aggregation mode 32 at the edges of the scan. The HSIs are 
projections of ASIs on the ground. It is sometimes referred to as ground sample distance (GSD) or ground sample 
interval (GSI) 191 . They vary with satellite altitude and scan angle. For SNPP, the variation of altitude is tightly 
controlled by drag make-up maneuvers to maintain 101.5 minutes orbital period and 16-day repeat cycle 171 . Over the 
equator, the altitude is within 1 km of 829.8 km. However, the altitude is about 840 km over the north pole and about 
855 over the south pole 57 ’ 1 . At the mean altitude of 840 km, the post-aggregation HSIs for M-bands are 786 m and 742 
m in the scan and track direction at nadir, respectively, and 1711 m and 1658 m at the edges of scan. For 1-bands, the 
HSIs are twice as finer as the HSIs for the M-bands. For the DNB at this mean altitude, HSIs are maintained within 7% 
of 770 m in the scan direction and within 3% of 765 m in the track direction throughout the scan. 

Figure 1 depicts four examples of scan direction LSFs. The LSFs for band 13 are typical for the 5 I-bands and the LSFs 
for band M6 are typical for the 16 M-bands. The M-band exception is band Mil which has a side-lobe about 5% the 
magnitude of the main lobe about 2.3 sampling intervals away in the non -aggregated zone. The dashed lines in the 
subplots for 13 and M6 represent the reference LSFs assuming the detector size and integration drag are equivalent to 
ASI without optical blur. LSFs for the DNB are mostly squares, with the effects of optical blur on the LSFs near the 
edges of scan. The two dashed lines represent LSFs in aggregation modes 18 and 31without optical blur. 
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Figure 1. Typical line spread functions (LSFs) for I-bands (represented by 13 here) and M-bands (represented by 
M6 here). LSFs in Agg2xl and Agg3xl zones are derived from LSFs in non-aggregated zones from ground testing. 

DNB LSFs are from HGA only for this paper. 


The LSFs in the scan direction may be parameterized in dynamic field of view (DFOV) and horizontal spatial 
resolution (HSR). DFOV is the full-width at half maximum (FWHM) of the LSF in the scan direction including the 
effect of time integration. It may be represented in term of ASI. The projection of DFOV on the ground is the ground 
DFOV (GDFOV). HSR is the half wavelength on the ground where modulation transfer function (MTF) of LSF drops 
to '/2 as expressed in Equation (3) below. It is sometimes regarded as “effective” detector (or sensor) resolution and may 
be expressed in terms of HSI. MTF is a normalized Fourier transfer function commonly used to quantify image quality, 
i.e., 


MTF {(D) = 


\FourierTansform{LSF)\ 

\FourierTansform{LSF)\ 


( 1 ) 


where CO is expressed in fraction of Nyquist frequency that is the maximal detectible sinusoidal waves by the Nyquist 
sampling theorem 1101 , i.e., 
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J Nyquist 

From the MTF curve, the FISR may be calculated as 


' 1 ' 

v 2 HSI j 


MTF 


r 1 " 

2 HSR 


= 0.5. 


From the LSF, we may also calculate the ensquared energy 111 , 


( 2 ) 


(3) 


EE = j* l LSF(x)dx (4) 

1 1 2 

which represents the total weight within the ASI. Integrated out-of-pixel (IOOP) spatial response may also be 
introduced to represent the total weight contaminated by the neighboring pixels, i.e., 


IOOP = f , LSF(x)dx = 1 -EE. (5) 

Jl x\> 

I 2 

To simplify these calculations of the scan direction EE and IOOP, we assume that the LSFs in the track direction are 
uniformly distributed within the sampled pixel (i.e., track EE = 1 and IOOP = 0). 

For example, an ideal triangular LSF with a FWHM the same as its ASI, the HSR is 1.13 HSI, MTF is 0.40 at the 
Nyquist frequency, EE is 75% and IOOP is 25%. For an ideal square LSF with a FWHM the same as its ASI, the HSR 
is 0.83 HSI, MTF is 0.63 at the Nyquist frequency, EE is 100% and IOOP is 0%. 

2.1 Sensor (detector) Resolutions 

Figure 2 shows the scan direction GDFOV to HSI ratio for each detector for the I-bands and M-bands in the three 
aggregation zones, as well as for the three DNB (HGA) aggregation modes (31, 22 and 18). For all 5 I-bands, pixels are 
under sampled at the top of atmosphere. The under-sampling is most acute for bands 13 to 15 in the Agglxl zones and 
may cause aliasing 1 ^ 1 . However, when they measure radiance from the ground surface, atmospheric modulation^ 121 may 
reduce the under-sampling effect as the optical path lengthens in the Agglxl zones, which are the furthest from nadir. 
For M-bands, pixels in the Agglxl zones are over-sampled. However, the pixels in the Agg2xl and Agg3xl zones are 
sampled at the Nyquist frequency. This Nyquist frequency sampling is also true for DNB aggregation modes 31, 22, 
and 18, and may be extrapolated to all DNB aggregation modes throughout the scan. 

Figure 3 shows the scan direction HSR to HSI ratio for each detector for the I-bands and M-bands in the three 
aggregation zones, and the three DNB (HGA) aggregation modes. By comparing Figures 2 and 3, we can see that for I- 
bands and M-bands in the Agglxl zones, the sensor HSR/HSI ratio is larger than the GDFOV/HSI ratio. However, as 
the aggregation zones approach nadir with more samples aggregated, the LSF becomes more square-like and GDFOV 
approaches 1.0 HSI, but HSR approaches 0.83 HSI. This is also true for DNB. Thus, the HSR for a square-like LSF 
may mis-represent “effective” sensor resolution as the footprint on the ground. In general, GDFOV is a better parameter 
representing sensor footprint, even though long tails and side- lobes, such as the one for band Mil, may not be well 
represented by GDFOV. 
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1.4 



Figure 2. Scan direction GDFOV to HSI ratio for each detector for the I-bands, M-bands and DNB (FIGA). 



Figure 3. Scan direction HSR to HSI ratio for the each detector for the I-bands, M-bands and DNB (HGA). 


2.2 Imaging Quality Parameters 

Image quality is commonly measured by modulation transfer function (MTF), a measure of contrast reduction of 
sinusoidal functions. Figure 4 shows the scan direction MTF for each detector for the I-bands and M-bands in the three 
aggregation zones at 1.0 and 0.5 Nyqusit frequencies. It also shows DNB (HGA) in three aggregation modes, which has 
values comparable to those for bands II and 12. The MTF for M-bands has lower values, especially for the visible bands 
Ml to M6. The effect of side-lobe in Mil LSF does not show up in the MTF at the 1.0 Nyquist frequency in the non- 
aggregation zones but does so at the 0.5 Nyquist frequencies. This demonstrates the need to quantify MTF at multiple 
frequencies. 

Figure 5 shows scan direction ensquared energy (EE) for each detector for the I-bands and M-bands in the three 
aggregation zones, and three DNB (HGA) aggregation modes. I-bands collect more energy within the sampled pixels 
than M-bands. As aggregation zone approaches nadir, more energy is collected within the sampled pixel. It may be 
inferred that DNB collects 100% energy within the sampled pixels as the aggregation modes approaches nadir 
(aggregation mode 1). 
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Figure 4. Scan direction MTF for the each detector for the I-bands, M-bands and DNB (HGA) at the 1.0 and 0.5 

Nyquist frequencies. 
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Figure 5. Scan direction Ensquared energy (EE) for the each detector for the I-bands, M-bands and DNB (HGA). 

Figure 6 shows the scan direction integrated out-of-pixel (IOOP) spatial response for each detector for the I-bands and 
M-bands in the three aggregation zones, and three DNB (HGA) aggregation modes. This is the amount of energy 
contaminated by the neighboring pixels. As we can see, band Ml pixels are contaminated the most, up to 35% from the 
neighboring pixels in the Agglxl zones. 
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Figure 6. Scan direction Integrated Out-Of-Pixel (IOOP) spatial response for the each detector for the I-bands, M- 

bands and DNB (HGA). 

2.3 On-orbit spatial response characterization 

For on-orbit VIIRS bands spatial response characterization (for verification purposes), we devise an approach that takes 
advantage of an available high contrast linear feature to provide a direct, although approximate, measurement of the 
VIIRS instrument along-scan LSF. In our approach we use a co-geo located Landsat TM image to precisely locate the 
linear feature in the VIIRS image. We use the Lake Ponchartrain Causeway as a high contrast linear feature. Figures 
7(a) and 7(b) show sections of Landsat TM and VIIRS images co-geolocated to a spatial resolution 13 times finer than 
the VIIRS image resolution (and slightly finer than the Landsat TM resolution). The data dropouts in the Landsat TM 
image are due to the well-known scan line corrector failure on Landsat 7 [13] (and are not included in the analysis). We 
locate the causeway in the Landsat TM image by processing the co-geolocated Landsat TM image with the HSeg [14] 
image segmentation and selecting region objects corresponding to the causeway with HSegViewer utility. We then 
compute the linear best fit equation to the found causeway pixel locations as a precise model of the causeway location. 
A plot of this best fit linear best fit equation is shown as the red line in Figure 7(c). We also noted where the center of 
each VIIRS image pixel was located in relation to the modeled causeway location, shown as the cyan colored dots in 
Figure 7(c). 



Figure 7. Co-geolocated Lake Ponchartrain subsets for (a) a Landsat TM image from 21 Nov. 2012 and (b) a VIIRS 
image from 18 Nov. 2012. (c) 16-times magnified plot of the linear equation fit to the northern end of the causeway 
(red line) and the locations of the centers of the SNPP VIIRS image data pixels (cyan dots). 

We then used the information embodied in Figure 7(c) to find the fractional column causeway location in each scan line 
of the VIIRS image that crossed the causeway, and plotted the VIIRS radiance values of the 5 pixels in each scan line 
closest to the causeway. For this example in the Agg2xl zone, we then fit these values to a sum of two identical 


Proc. of SPIE Vol. 8866 88661 G-7 


Downloaded From: http://proceedings.spiedigitallibrary.org/ on 10/12/2013 Terms of Use: http://spiedl.org/terms 




Gaussians curves (each centered at +/- 0.5 HSI). The LSF measurements, model fit and pre-launch measurement are 
shown in Figure 8 for bands 12 and 13, the raw data points and the modeled curves tightly follow the pre -launch curves. 
For bands II and 14 there is a much larger scatter in the data points. The most scattered data points can be found in 
band 15. However, the modeled curves are still similar to the pre-launch curves. 
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Position (post-aggregation HSI) Position (post-aggregation HSI) Position (post-aggregation HSI) 


Figure 8. Scan direction LSFs for 5 I-bands in Agg2xl zone, measured by the bridge, the double-Gaussian fit and 

the pre- launch data. 


These modeled LSFs yield FWHM as the sensor GDFOVs for the 5 I-bands. The mean values and standard deviations 
binned in the three aggregation zones are plotted in Figure 9, along with pre-launch data. As we can see, the on-orbit 
LSFs measurements based on the ground target compare favorably with the pre -launch test, especially in the 
aggregation zones close to the start of scan, especially considering the atmospheric modulation [12] . For the GDFOV 
measured from the end of the scan, the sensor zenith angles are in the opposite direction from the solar zenith angles 
and so sun glint may affect the GDFOV retrievals. Further study of the effects of sun glint is on-going, along with the 
retrievals of GDFOVs for the M-bands and DNB. 
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Figure 9. Scan direction GDFOV to HIS ratio for the 5 I-bands using Gaussian fit as shown in Figure 8 in the 5 
aggregation zones from the start to end of scan. 
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3. BAND-TO-BAND CO-REGISTRATION (BBR) 

The BBR is a statistical value of a collection of individual detector-to-detector co-registration (DDR) measurments 
between two corresponding spectral bands. DDR may be defined as the overlap of point spread functions (PSFs) from 
one detector in one spectral band to a corresponding detector in another spectral band [4,15 ’ 16 l 

DDR, j = 1 - g \PSF x - PSFj | / 2 dxdy (6) 

(x,y)f=R 


where PSF is normalized, i.e., 

g PSF dxdy = 1 . (7) 

(x,y)<ER 

Since the exact PSF is usually difficult to obtain precisely, especially for space-borne sensors which are subject to 
various conditions, an approximation of Equation (6) is used in practice. Measurements are done to estimate apparent 
relative locations (offsets) of detectors between corresponding band pairs. The offsets themselves are then be used as 
measures of BBR, as is done for the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard NASA Earth 
Observing System (EOS) Terra and Aqua satellites 117 ' 18 ' 19 ' 201 . For VIIRS, the offsets used to compute BBR matrices are 
based on a chosen reference area 14 with the assumption that an equivalent PSF is evenly distributed over that area. This 
may be expressed as 


DDR iJ «( (8) 

L s L T 

where AS, ■ and AT, . are the measured detector offsets (in HSI units) and L s and L, are the length scales 
approximating the equivalent area in the scan and track directions, respectively. 

The VIIRS performance requirements specify the length scales at “specified HSR” at the ends of scan in both scan and 
track directions at 1600 m for M-bands and at 800 for I-bands. As we described in the beginning of Section 2, such 
length scales are comparable with the HSI at the ends of the scan. Thus, the HSI was used as the length scales to assess 
whether VIIRS met its BBR requirement^ 1 . In the following subsections, we are going to present the offsets measured 
pre-launch to gauge the “improvements” of BBR in the Agg2xl and Agg3xl zones relative to the Agglxl zone BBR 
using the HSI as defined in the VIIRS requirements. We will then compare this to the BBR from alternative approaches 
based on the scan direction length scales defined by the GDFOV and HSR. In this analysis, the length scale in the track 
direction is kept the same as track HSI. Since this paper focuses on the comparisons purpose, we ignore the effects of 
on-orbit conditions (e.g. atmosphere blurring). 

3.1 Band-to-band offsets 

Figure 10 shows the band averaged offsets relative to band II in the scan and track directions. Data for 2x1 and 3x1 
aggregation zones are derived from measurements in thermal vacuum at the nominal temperature performance plateau 
in the diagnostic mode, i.e., Agglxl zone. The relative positions of all detectors were precisely measured 141 . We only 
show band averaged values here for clarity. 
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Figure 10. Band average offset from ground test under thermal vacuum condition at the nominal temperature 

performance plateau. 


3.2 BBR matrices based on HSI 

In the Agglxl zone the DDRs have been calculated based on the HSIs in the scan and track directions. Note, since 
these BBR results are for the Agglxl zone in which 4 (8) M-band (I-band) detectors are deleted on-board, these DDRs 
are based on the center 12 (24) M-band (I-band) detectors. Taking the minimum of the DDRs as the representation of 
BBR between the band pairs, the BBR matrix (in %) for every band pair is presented in the upper-right triangle in 
Figure 11. The lower-left triangle is the VIIRS BBR specification for “at least 99.7% of corresponding pixel samples”. 
To make the figures easier to understand, colors in the matrix indicate how well each band-pair is co-registered (with 
“green” the best and “pink” the worse), they are not meant to indicate compliance with the requirement. All VIIRS 
band-pairs met the requirement except pairs M9/M13 and M11/M13, which were only non-compliant by a small 
amount. 

We then use the offsets in term of post-aggregation HSI and compute the BBR matrices in the Agg2xl and Agg3xl 
zones, as shown respectively in the upper-right and lower-left triangles in Figure 12. In this figure, we see the 
“improvements” of BBR up to 10% and 13% respectively in the Agg2xl and Agg3xl aggregation zones. 
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|< - BBR Sepc for 99.7% of corresponding pixel samples based on "specified HSR 800 m or 1600 m end of scan > 

70 : Spec 97 :>94% 91 : > 89% 87 :>84% 82 :>79% 78 :<80% 


Figure 1 1 BBR based on HSI. Lower left triangle: specification at 99.7% confidence level. Upper right triangle: 
minimum among detector pairs based on HSI in the Agglxl zones. 
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|< BBR = Min( 16/32 HSI-based DDRs for HAM sides A&B) in 3x1 aggregation zones— — >| 


97 : > 94% 91 : > 89% 87 : > 84% 

Figure 12 BBR based on the HIS in Agg2xl zones (upper right triangle) and Agg3xl zones (lower left triangle). 
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3.3 Alternative parameterization of BBR based on GDFOV and HSR 

As we shown in Figure 2 and Figure 3, the GDFOV and FiSR in the scan direction are different from the HSI, 
especially in the Agglxl zones. As above, in this section we set the length scale to the HSI (which is about the same as 
GIFOV [4 ^), but change the scan direction length scale to GDFOV and HSR. The results are shown in Figure 13. 
Compared to the BBR matrix based on HSI, the BBR based on GDFOV “improves” (up to 4%) for the M-band pairs 
but “degrades” (up to 3%) for the I-band pairs, while BBR based on scan HSR “improves” the M-band pairs (up to 4%). 
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Figure 13 BBR matrices based on scan GDFOV (upper right) and scan HSR (lower left). 

To be consistent, when we use the scan HSR for the length scale in the scan direction, we should use track HSR for the 
length scale (which is close to the theoretical value of 0.83 GIFOV^ 4 ') for the track direction as well. This is done and 
the resultant BBR matrix is compared with the BBR matrix based on HSI (Figure 14). Using HSR as the base for both 
the scan and track direction, the BBR “improves” (up to 4%) for the M-band pairs but “degrades” (up to 3%) for the I- 
band pairs. Using HSR as the length scale for nearly square LSFs, the offset between the band pairs are amplified 20%. 
For the case of SNPP, the offsets in the track direction are relatively small (Figure 10) and thus the amplification is 
small. If the offsets are relatively large (similar to scan offsets in Agglxl zone as shown in Figure 10), then, that 20% 
amplification could eat away the small BBR margin. 
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Figure 14 BBR matrices based on HSI (upper right) and scan/track HSR (lower left). 

In theory, for an ideal triangular LSF (which is a close approximation of scan LSF [5 ' 17 ' 19 ' 20 ') with an FWHM the same as 
its ASI and GDFOV the same as its HSI, the HSR is 1.13 HSI, 13% larger than the footprint usually represented by the 
FWHM. For an ideal square LSF (which is a close approximation of track LgpHs, 18,19,20]^ an FWHVI the same as 
its ASI and GIFOV the same as its HSI, the HSR is 0.83 HSI, 17% smaller than the footprint usually represented by the 
FWHM. HSR is calculated as the half wave-length where modulation transfer function (MTF) of LSF drops to 'A, i.e., 
the contrast of a sinusoidal wave at that wave-length is reduced to 'A. This is sometime called effective resolution. It 
lacks the general meaning of representing the size of the PSF or LSF patterns projected on the ground. It could be mis- 
leading to use HSR as a measure of image quality or as a length scale to calculate BBR. Therefore, HSR as defined in 
Equation (3) should not be used in the field of remote sensing, especially when the remote sensed data is used to further 
retrieve biogeophysical parameters. 


4. CONCLUDING REMARKS 

In this paper, we present the parameterization of sensor spatial responses and band-to-band co-registration (BBR) in 
various ways. These results are derived from recently re -processed data measured before the SNPP VIIRS was 
launched in October 2011. For the VIIRS spatial response in the scan direction, all 5 I-bands are under-sampled as 
measured by their ground dynamic field of view (GDFOV). The image quality as measured by modulation transfer 
function (MTF), ensquared energy (EE) and integrated out-of-pixel (IOOP) spatial response are high. However, under- 
sampling may cause aliasing. For the 16 M-bands, GDFOV values in the scan direction indicated that pixels in the non- 
aggregated zones near the edges of scan are over-sampled, their image quality measured by MTF is relatively low 
(compared to I-bands) and contamination from neighboring pixels as indicated by IOOP is relatively high (compared to 
I-bands). In the Agg2xl and Agg3xl zones, M-bands pixels are sampled to the Nyquist frequency. In terms of BBR, all 
band pairs have offsets less than 20% of their corresponding pixel sizes. They have near 80% or better footprint 
overlaps, especially in the Agg2xl and Agg3xl zones. When comparing HSI, GDFOV and HSR as the length scales for 
fractional band-to-band mis-registration, we found that GDFOV calculated as the full-width at half-maximum (FWHM) 
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of sensor spatial response function is a better representation of footprint size. We also propose that the HSR as defined 
in Equation (3) is neither an adequate representation of the size of sensor spatial response nor an adequate measure of 
imaging quality, and therefore HSR concept should not be used in the field of remote sensing when the remote sensed 
data is used to retrieve biogeophysical parameters. 
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